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Abstract 



We present a full diagrammatic computation of the one-loop corrections from the neu- 
trino/sneutrino sector to the renormalized neutral CV-even Higgs boson self-energies 
and the lightest Higgs boson mass, M^, within the context of the so-called MSSM- 
seesaw scenario. This consists of the Minimal Super symmetric Standard Model with 
the addition of massive right handed Majorana neutrinos and their supersymmetric 
l/"~) ■ partners, and where the seesaw mechanism is used for the lightest neutrino mass gen- 

ial . eration. We explore the dependence on all the parameters involved, with particular 

emphasis in the role played by the heavy Majorana scale. We restrict ourselves to the 
case of one generation of neutrinos/sneutrinos. For the numerical part of the study, we 
consider a very wide range of values for all the parameters involved. We find sizeable 
corrections to Mh, which are negative in the region where the Majorana scale is large 
(10 13 — 10 15 GeV) and the lightest neutrino mass is within a range inspired by data 
(0.1 — 1 eV). For some regions of the MSSM-seesaw parameter space, the corrections 
to Mh are substantially larger than the anticipated Large Hadron Collider precision. 
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1 Introduction 

The current impressive experimental data on neutrino mass differences and neutrino mixing 
angles [1] indicate clearly a signal of new physics beyond the so far successful Standard 
Model of Particle Physics (SM). In order to incorporate the non- vanishing neutrino masses 
required by data an extension of the SM with massive neutrinos is mandatory. Among the 
various possibilities to extend the SM we choose here the most popular one that incorporates 
massive Majorana neutrinos and also stabilizes the electroweak symmetry breaking scale, 
v = 174 GeV, against potentially large radiative corrections in the presence of the new 
physics scale. We refer to the simplest version of a supersymmetric extension of the SM, 
the Minimal Supersymmetric Standard Model (MSSM) |2], with the addition of heavy right- 
handed Majorana neutrinos, and where the well known seesaw mechanism of type I [3] is 
implemented to generate the observed small neutrino masses. From now on we will denote 
this model by "MSSM-seesaw" . 

In this MSSM-seesaw context, the smallness of the light neutrino masses, m„ ~ m D /mM, 
appears naturally due to the induced large suppression by the ratio of the two very distant 
mass scales. Namely, the Majorana neutrino mass tum, that represents the new physics scale, 
and the Dirac neutrino mass mp, which is related to the electroweak scale via the neutrino 
Yukawa couplings Y v , by mo = Y v v sin/3. The Higgs sector content in the MSSM-seesaw 
is as in the MSSM [4], and tan/3 is given, as usual, by the ratio of the two MSSM Higgs 
vacuum expectation values (v.e.v.s). Although the present neutrino data requires two or 
more neutrino generations, we shall adopt here the simplest case of one neutrino generation 
in order to fully understand first the role of one single Majorana scale rriM, and postpone the 
more complex case of three generations for a future work. In this simplified one-generation 
MSSM-seesaw framework, small neutrino masses of the order of m v ~ 0.1 — 1 eV can be 
easily accommodated with large Yukawa couplings, Y v ~ 0(1), if the new physics scale is 
very large, within the range tum ~ 10 13 — 10 15 GeV. This is to be compared with the Dirac 
neutrino case where, in order to get similar small neutrino masses, extremely tiny, hence 
irrelevant, Yukawa couplings of the order of Y u ~ 10 -12 — 10~ 13 are required. 

The hypothesis of Majorana massive neutrinos is very appealing for various reasons, 
including the interesting possibility of generating satisfactorily baryogenesis via leptogene- 
sis |5], and also because they can produce an interesting and singular phenomenology due 
to their potentially large Yukawa couplings to the Higgs sector of the theory, the MSSM in 
the present case. Among the most striking phenomenological implications of these MSSM- 
seesaw scenarios [5], it is worth mentioning: 1) the prediction of sizeable rates for lepton 
flavor violating processes, indeed within the present experimental reach for specific areas of 
the model parameters [3(8], 2) non-negligible contributions to electric dipole moments of 
charged leptons [9], and 3) the occurrence of sneutrino-antisneutrino oscillations [TU] and 
sneutrino flavor-oscillations [TTj . 

The present paper investigates another implication of heavy Majorana neutrinos that 
could be as relevant as these previously mentioned ones. More specifically, we are interested 
here in the indirect effects of Majorana neutrinos via their radiative corrections to the MSSM 
Higgs boson masses. In particular, our study will be focused on the radiative corrections 
to the lightest MSSM CP-even h boson mass, M^, due to the one-loop contributions from 
the neutrino/sneutrino sector within the MSSM-seesaw framework. Previous studies in par- 



ticular SUSY scenarios and under specific assumptions on the model parameters J8J 
indicate that the size of these radiative corrections to the Higgs mass parameters in the case 
of extremely heavy Majorana neutrinos can be sizeable due to the large size of Y v . 

For the estimates of the total corrections to M^ in the MSSM-seesaw, obviously, the one- 
loop corrections from the neutrino/sneutrino sector that we are interested here have to be 
added to the existing MSSM corrections. The status of radiative corrections to M^ in the non- 
vjv sector, i.e. in the MSSM without massive neutrinos, can be summarized as follows. Full 
one-loop calculations [15] have been supplemented by the leading and subleading two-loop 
corrections, see [16] and references therein. Together with leading three-loop corrections [17] 
the current precision in M^ is estimated to be ~ 2 — 3 GeV [TB] . 

Regarding the previous estimates of neutrino/sneutrino radiative corrections to the Higgs 
mass parameters the status is as follows. In Ref. [12] the one-loop corrections to Mh were 
estimated within a split SUSY scenario where the soft-SUSY-breaking mass associated to 
the right handed neutrino, m^, was chosen to be very large, of the order of the Majorana 
scale rriM- They worked in the zero external momentum approximation and switching off the 
SU(2) x U(l) gauge interactions. Besides, they used the mass insertion approximation for the 
other soft-breaking sneutrino parameters, A v and B u , associated to the trilinear coupling and 
neutrino 5-term respectively. A large and negative correction from the neutrino/sneutrino 
sector of the order of a few tens of GeV was found for rriM = 10 14 GeV and m^ ~ 0{mu)- In 
Ref. [13] the radiative one-loop effects of the neutrino S-term on the Higgs mass parameters 
within the context of mSUGRA (with universal scalar masses at the mgusY? including m^) 
were analyzed by means of the renormalization group equations (RGEs). They found large 
effects from this B v term that indeed could destabilize the electroweak symmetry breaking. 
By requiring a proper breaking in this mSUGRA framework they concluded with an upper 
bound of B U Y 2 / (8n 2 ) < msusy/tan/3. Large corrections to the Higgs soft mass parameters 
within a SUSY-seesaw framework with total or partial universality conditions have also been 
found by a similar RGEs analysis in [8j[Il]. In [8] it was concluded that these corrections 
induce a considerable decrease in the physical Higgs boson masses which in turn enhance the 
rates of the Higgs- mediated LFV processes. In [2] the large threshold corrections found from 
the heavy neutrinos/sneutrinos were shown to affect, and even dominate at large B u , the 
radiative breaking of the electroweak symmetry and also modify considerably the predictions 
on the neutralino dark matter abundance. 

In this work, we will consider instead the more general MSSM-seesaw scenarios with no 
universality conditions imposed, and explore the full parameter space, without restricting 
ourselves just to large or small values on neither of the relevant neutrino/sneutrino param- 
eters. In principle, since the right handed Majorana neutrinos and their SUSY partners are 
SU(2) x U(l) singlets, there is no a priori reason why the size of their associated parameters 
should be related to the size of the other sector parameters. In the numerical estimates, we 
will therefore explore a wide interval for all the involved neutrino/sneutrino relevant input 
parameters. 

We will present here a full one-loop computation of the radiative corrections to the lightest 
CP-even Higgs boson mass from the (one generation) neutrino/sneutrino sector in which we 
will not use any of the previous approximations and we will not set the external momentum 
to zero. The complete set of one-loop neutrino/sneutrino contributing diagrams will be taken 
into account, with both Yukawa and gauge couplings switched on. We also analyze the results 



in several renormalization schemes, which will be shown to provide remarkable differences. 
In addition, we present some analytical and numerical results in the interesting limit of very 
large mu as compared to all other scales involved, which will help us in the understanding 
of the important issue of the decoupling/non-decoupling of the heavy Majorana scale. Our 
further study in the particular region of large tum and m^ will also allow us to compare our 
results with those in [12] . 

Our final aim is to find out to what extent the radiative corrections computed here enter 
into the measurable range. The experimental perspectives for the Higgs mass measurements 
with precision enough to be sensitive to such sizeable radiative corrections, as the ones 
found here, are indeed quite promising. The LHC has good prospects to discover at least 
one neutral Higgs boson over the full MSSM parameter space and a precision on the mass of 
a Standard Model (SM)-like Higgs boson of ~ 200 MeV are expected pH2T] (see e.g. [22l|23] 
for reviews). At the ILC a determination of the Higgs boson properties (within the kinematic 
reach) will be possible, and an accuracy on the mass could reach the 50 MeV level [24Tl2~T] . 
The interplay of the LHC and the ILC in the neutral MSSM Higgs sector will improve 
certainly these measurements [28|[29]. 

The paper is organized as follows. In section 2, we summarize the most important 
ingredients of the MSSM-seesaw scenario that are needed for the present computation of 
the Higgs mass loop corrections. These include, the setting of the model parameters and 
the complete list of the Lagrangian relevant terms. A complete set of the corresponding 
relevant Feynman rules in the physical basis is also provided here. They are collected in the 
Appendix A and, to our knowledge, they are not available in the previous literature. We 
also comment shortly in section 2 on the comparison between the Dirac and the Majorana 
cases. In section 3 we present the renormalization procedure and emphasize the differences 
between the selected renormalization schemes, specifically, the on-shell and the DR schemes. 
Section 4 is devoted to the results. First we present the analytical results for the renormalized 
Higgs boson self-energies (the main formulas are collected in Appendix B). Then we present 
the numerical results in terms of all the relevant neutrino/sneutrino parameters that we 
explore exhaustively in the full plausible range. We also include in this section a study of 
the behavior of the renormalized Higgs self-energies in the large mu limit. The final part 
of this section summarizes the main numerical results for the lightest Higgs boson mass 
corrections. Finally, section 5 contains the conclusions. 

2 The MSSM-seesaw model 

The model we are interested in here is the MSSM extended by right handed neutrinos and 
their SUSY partners, and where a seesaw mechanism of type I [3] is implemented to generate 
the neutrino masses and mixing angles. This is called usually the MSSM-seesaw model. For 
simplicity, as already announced in the introduction, we will restrict here to the one genera- 
tion neutrinos/sneutrinos case although the full compatibility with present neutrino data for 
mass differences and mixing angles, requires additional neutrino generations. Since the main 
idea is to analyze the radiative corrections from the neutrino-sneutrino sector to the lightest 
Higgs mass, we restrict ourselves to the case of one generation of neutrinos/sneutrinos. We 
illustrate first this simpler case and postpone the more complex case of three generations for 



a future work. 



2.1 The neutrino/sneutrino sector 

The MSSM-seesaw model with one neutrino/sneutrino generation is described in terms of 
the well known MSSM superpotential plus the new relevant terms contained in: 



W = e l3 



Y V R\UN - YiH[L j R 



l -Nm M N, (1) 



where tum is the Majorana mass and N = (D R , (vr) c ) is the additional superfield that 
contains the right-handed neutrino v R and its scalar partner v R . Here and in the following f c 
denotes the particle-antiparticle conjugate (c-conjugate in short) of a fermion / (f c = Cf T ) 
and /* denotes the complex conjugate of sfermion /. The lepton Yukawa couplings are Yi >v , 
and we use the convention €12 = — 1 • The other superfields, L containing the lepton (ul, e£) 
and slepton (ul,gl) 577(2) doublets, R containing the lepton (eR) c and slepton e* R SU(2) 
singlets, and H\p containing the Higgs boson 577(2) doublets and their SUSY partners, are 
as in the MSSM. We follow here the notation of [I]. 

There are also new relevant terms in the soft SUSY breaking potential due to the addi- 
tional sneutrinos vr [TO] : 



Koft = m \ v l v L + m \vRVR + {Y v A v H\v L v R + m M B v VRVR + h.c.) . (2) 

After electro-weak (EW) symmetry breaking, the charged lepton and Dirac neutrino 
masses can be written as 

m\ = Yi v 1} m D = Y u v 2 , (3) 

where Vi are the vacuum expectation values (VEVs) of the neutral Higgs scalars, with vu2) = 
v cos(sin)/3 and v = 174 GeV. 

The 2x2 neutrino mass matrix is given in terms of mjj and rriM by: 

W = ( ° mD V (4) 

^ m D m M ) 

Diagonalization of M v leads to two mass eigenstates, rti (i = 1,2), which are Majorana 
fermions: 






m = v — cos6(u L + (isl) c ) - sin.0(vR + {vr) 

n 2 = N = sin 9(y L + (u L ) c ) + cos 9{v R + {v R ) c ) (5) 



with the respective mass eigenvalues given by: 

1 



™>v,N=2\ mM ^ V m M + 4m D J • ( 6 ) 

It should be noticed that we have introduced an alternative notation that makes it easier to 
identify the specific neutrino by its mass: v is the lighter one and N is the heavier one. It 



should also be kept in mind that with this convention m v < and m^ > 0, but the physical 
Majorana neutrino states have the proper positive masses. These physical neutrinos can 
be reached by an additional rotation, v — > z/ = e* 757r / 2 z/ = —i^^u, leading to m u i = \m v \. 
However, we prefer to work instead with the mass eigenstates in (J3J) to avoid extra % and 
75 factors in the computation. Of course the final results in this work for the Higgs mass 
corrections are not sensitive to this choice. 

The mixing angle that defines the mass eigenstates is given by, 

tan£ = -^ = ^. (7) 

m D m N 

Other useful relations between the model parameters mp, txim and the physical neutrino 
parameters, m„, m^ and 9 are the following: 



sin 2 # 




m M 



'mr M + 4m 2 D 



2 a m N l L , m M \ ,„s 

cos 9 = = - 1 H = , (9) 

m N -m v 2 \ y/mlf + Am 2 D J 

m D = -\f{m N - m v ) 2 - (m N + m u ) 2 , (10) 

ml = -m u m N , (11) 

m M = m u + m N . (12) 

Regarding the sneutrino sector, the sneutrino mass matrices for the CP-even, M + , and 
the CP-odd, M_, subsectors are given respectively by [TU] : 

X>2 __ ( m 2 L + m 2 D + |M|cos2/3 m D {A v - fi cot (3 ± m M ) \ ,^ 

± \ m D (A u - jj, cot f3 ± 771m) fn\ + rn 2 D + m 2 M ± 2B v m M ) 

The diagonalization of these two matrices, M|, leads to four sneutrino mass eigenstates, 
fti (i = 1, 2, 3, 4) with respective CV parities CV{ri\^) = +1 and CP(n 3i4 ) = — 1: 

ni = 1/+ = v2(cos 9 + Re z>^ — sin 9 + Re z>r) , 
n 2 = N + = V2(sin 9 + Re u L + cos # + Re i> R ) , 

77.3 = z>_ = v2(cos 9- Im ^ — sin 9_ Im Pr) , 

77.4 = iV_ = v2(sin 6L Im Vl + cos #_ Im z/r) . (14) 

It should again be noted that we have introduced an alternative notation that makes it easier 
to identify the specific sneutrino by its parity and mass: z> + , iV + are respectively the lighter 
and the heavier ones with CV = +1, and z>_, N_ are the lighter and the heavier ones with 
CV — — 1. The corresponding mass eigenvalues are: 

m l+,N + = 2^ + m ! + m | + 2m l + 2 M z cos 2/3 + 2B - m M ) (15) 

T g J 4m D( A » - Ai cot /3 + m M ) 2 + (m 2 f + m\ - m\ - -M 2 cos 2/3 + 2B u m M ) 2 , 
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5_,iV_ 



-(m 2 M + m\ + m\ + 2m 2 D + -M§ cos 2/3 - 2B u m M ) 



(16) 



T o\/ 4m i)(^ -/icot/3-m M ) 2 



m 



,-' M + m| - m| - -Mf cos 2/3 - 2B v m M ) 

The mixing angles in the two subsectors are given respectively by: 

2m£ ) (A l/ — n cot /3 ± m M ) 



sin 2# H 



Am 2 D (A v — /icot/3 ± m A /) 2 + (m| f + m| — m~ — \M\ cos 2/3 ± 2B u itim) 



;i7) 



2.2 The Higgs boson sector at tree-level 

In this subsection we summarize the Higgs-boson sector of our model at tree-level. Contrary 
to the SM, in the MSSM two Higgs doublets are required. The Higgs potential [30] 



V = m\\Ui\ 2 + ml\U 2 \ 2 -ml 2 (e ab UlU b 2 + h.c.) 

+ \{9 2 + 9 2 ) [l^il 2 - l^l 2 ] 2 + \9 2 HH 2 \ 2 , (18) 

contains mi,m 2 ,mi 2 as soft SUSY breaking parameters; g,g' are the SU(2) and U(l) gauge 
couplings, and e± 2 = —1. 

The doublet fields Hi and H 2 are decomposed in the following way: 



n 2 



ni J 


V -0r 




( <t>t 

\v 2 + j= 2 {4>l + ixl 



(19) 



The potential ( Tl8l) can be described with the help of two independent parameters (besides g 
and g')\ tan/3 = v 2 jv\ and M\ = — m^ 2 (tan/3 + cot /3), where M^ is the mass of the CP-odd 
Higgs boson A. 

The diagonalization of the bilinear part of the Higgs potential, i.e. of the Higgs mass 
matrices, is performed via the orthogonal transformations 



H 
h 
G 
A 
G ± 



cos a sin a 

— sin a cos a 
cos /3 sin f3 

— sin /3 cos /3 
cos /3 sin /3 

— sin /3 cos j3 



0? 

X? 

xl 



(20) 
(21) 
(22) 



The mixing angle a is determined through 

-(M| + M|)sin/3cos/3 



a = arctan 



M 2 cos 2 + Ml sin 2 /3 - m 2 



7T 



< a < 



(23) 



One gets the following Higgs spectrum: 



2 neutral bosons, CV = +1 

1 neutral boson, CV = — 1 

2 charged bosons 

3 unphysical Goldstone bosons 



h,H 

A 

H + ,H- 

Lr, Lt , Lt 



(24) 



At tree level the mass matrix of the neutral CP-even Higgs bosons is given in the 0i-02- 
basis in terms of Mz, Ma, and tan/3 by 



M 2 

1V1 Higgs 






M^ sin 2 (3 + M 2 Z cos 2 (5 - [M\ + M 2 ) sin (5 cos /3 
- (M 2 + M|) sin cos /3 M| cos 2 p + M\ sin 2 /3 

which by diagonalization according to fl20l) yields the tree-level Higgs boson masses 



(25) 



m 



H,h 



M\ + M 2 Z ± J (M 2 + Ml) 2 - AM 2 Z M 2 A cos 2 2/3 



The charged Higgs boson mass is given by 

m 2 H± = M 2 A + Mlf . 
The masses of the gauge bosons are given in analogy to the SM: 



(26) 



(27) 



M 



2/ 2 i 2\ 



m: 



-^(/ + ^)K + ^ 2 ); 



M 7 = 0. 



(28) 



2.3 The interaction Lagrangian 

Finally the interaction Lagrangian that is relevant for the present work, expressed in the 
{ v l, vr), [y L , v R ) electroweak interaction basis, is given by: 



r - r 



H 



c 



uZ + £-i>H 



c, 



(29) 



Here L v n and Lyu contain the interactions of the neutrinos and sneutrinos with the Higgs 
bosons respectively; and L v z and C^z those of the neutrinos and sneutrinos with the Z boson 
respectively 

For the various terms in (12"9"1) we find the following expressions: 



C, 



H 



-X( u lVr + v r i/ l )(H sin a + /i cos a) -i{v L v R - i/ R i/ L )AcosP) , (30) 



2M W sin j3 



4z 

C>vH 



- [(VZ>fu L )Zp] , 



2 cos 9w 



(31) 



gm D 



2M W sin p 
M w sin P 



H \iy*L v R + u lv r )(— H cosa + hsma)} 
[{v R v R + v* l vl)(H sina + h cos a)] 



. igm D _ 

[(i>* L i> L )(Hcos(a + (3) -hsin(a + /?))] 



2cos9w 
gm D 

2M W sin 



A u [(v* L i> R + vl v *r) (H sin a + hcosa)] 



gm D m M u ~ ~ . ~*~* V rr , . m 

' 2M w sM [( ^ + ^ R )(Hsma + hcosa)} 

. gm D m M r ,_ - ~*~ M , fl i 



^m£ 



4M 2 , sin 2 £ 
<? 2 



[(z/2z/ L ) (if 2 sin 2 a + /i 2 cos 2 a + A 2 cos 2 f3 + hH sin 2a 



8 cos 2 #vk 
g 2 m 2 D 



4M 2 , sin" /3 



[(z^l) (# 2 cos 2a - /i 2 cos 2a - A 2 cos 2/3 - 2hH sin 2a)] 

[(z>^ R )(F 2 sin 2 a + /i 2 cos 2 a + A 2 cos 2 /3 + /i# sin 2a)] , (32) 



'PZ 



1 9 



2 cos 9\v 



.■h$. 



,2 



|+ 4co ^ 2 ^ [(P2^)(^^)] • (33) 



The corresponding Feynman rules, expressed in the mass eigenstate basis, are collected in 
the Appendix A. Notice that this complete set of Feynman rules is, to our knowledge, not 
available in the literature so far. 

Some comments are in order. In the previous interaction Lagrangian, and consequently in 
the Feynman rules, there are terms already present in the MSSM. These are the pure gauge 
interactions between the left-handed neutrinos and the Z boson, given in (I3TJ) . those between 
the 'left-handed' sneutrinos and the Higgs bosons, given in fl32|) . and those between the 'left- 
handed' sneutrinos and the Z bosons, given in (133]) . In addition, in this MSSM-seesaw 
scenario, there are interactions driven by the neutrino Yukawa couplings (or equivalently 
m D since Y v = (gm D )/(\^2Mwsm(3)), and new interactions due to the Majorana nature 
driven by tum- These genuine Majorana terms are those in the seventh and eight lines of 
( )32|) and are not present in the case of Dirac fermions. 

2.4 Parameters and limits 

Regarding the size of the new parameters that have been introduced in this model, in addition 
to those of the MSSM, i.e., rriM, mo, m^, A v and B v , there are no significant constraints. In 
the literature it is often assumed that vim has a very large value, vn,u ~ (9(10 14-15 ) GeV, in 
order to get small physical neutrino masses \m u \ ~ 0.1 - 1 eV with large Yukawa couplings 
Y v ~ 0(1). This is an interesting possibility since it can lead to important phenomenological 
implications due to the large size of the radiative corrections driven by these large Yukawa 
couplings. In this paper we will explore, however, not only these extreme values but the full 
range for tum from the electroweak scale ~ 10 2 GeV up to ~ 10 15 GeV. 

On the other hand, the new soft SUSY-breaking parameters introduced in the sneutrino 



sector could be unrelated to those of the MSSM, or could be related, for instance, in the case 
one imposes (by hand) some kind of universality conditions. Whereas the non-singlet soft 
mass parameter mi, being common to the charged 'left handed' slepton, is constrained by 
the solution to the hierarchy problem to lie below a few TeV, the singlet soft mass m R is not, 
because it is not connected to the electroweak symmetry breaking at tree level. The other 
sneutrino soft mass parameters, B v and A v are not connected either. However, they can 
generate a mass-splitting between sneutrinos and antisneutrinos which in turn and via loop 
corrections can generate neutrino mass splittings [H] that are experimentally constrained. 
Then, if msusY represents a generic low SUSY breaking scale, with msusY ~ O(10 3 ) GeV 
one expects that \A V \, \B„\ < m SUSY [13j. According to these constraints, we will explore 
in this work values of these soft parameters ranging from the electroweak scale up to a few 
TeV. Besides, and due to the peculiarity of the behavior with m R and B v , as will be shown 
later, we will explore in addition the less conservative but interesting possibility where m R 
or B v are close to tum- 

For illustrative purposes and a clear understanding of our full one-loop results, three 
interesting limiting cases will also be considered in this work. 

(1) The seesaw limit: 

This assumes a large separation between the two neutrino mass scales involved, the 
Majorana mass and the Dirac mass, m^i 3> thd- Notice that both masses are different 
from zero, tjim 7^ and mp ^ 0, in this seesaw limit and, as we have said above, Y u can 
be large. The predictions are then given in power series of a dimensionless parameter 
defined as, 

e^«i. P4) 

m M 
The light and heavy neutrino masses are given in this limit by: 

2 
77? 

m v = -m D £ + 0(m D £ 3 ) ~ ^ , (35) 

m M 

m N = m M + 0(m D £) ~ m M . 

Furthermore, the mixing angle 9 is small in this limit and, therefore, v is made pre- 
dominantly of vl and its c-conjugate, {vl) c , whereas N is made predominantly of u R 
and its c-conjugate, (vr) c - 

In the sneutrino sector several mass scales are involved. Consequently, one has to set 
as an extra input their relative size to tum- The simplest assumption is to set the 
value of tum to be much larger than all the other mass scales involved, i.e., mu 3> 
rriD, Mz, /i, mi, m R , B v , A v . In this limit the sneutrino masses are given by: 

™l + ,u- = rn\ + -M 2 Z cos 2/3 t 2m D {A v - fj, cot f3 - B„)£ , 
m% ^_ = m 2 M ± 2B v m M + m\ + 2m 2 D . (36) 

The mixing angles 9± are small in this limit and, therefore, z> + and z/_ are made pre- 
dominantly of vl and its c-conjugate, v* L , whereas N + and AL are made predominantly 
of ur and its c-conjugate, v R . 



(2) The Dirac limit: 

In this limit one sets thm — (and tjid 7^ 0) and one recovers the neutrinos as any 
other fermion of the MSSM, i.e., as Dirac fermions. In the basis that we have used 
in ([5]) this is manifested by the fact that when rriM = 0, the two Majorana neutrinos v 
and TV are degenerate with m u = —thd and m^ = +m£>, and they combine maximally, 
i.e. with 9 = tt/4, to form a four component Dirac neutrino with mass m^. On the 
other hand, the sneutrino sector in this Dirac limit simplifies as well. When tum = 0, 
the real scalar fields get degenerate in pairs, 



m l+ = ml_ 



-(m| + m| + 2m^ + - Mf cos 2/3) (37) 



-\Lm 2 D (A v - /icot/3) 2 + (m 2 - - m 2 - - -Mf cos 2/3)' 

/ \l IX J-/ ) 



— (ml + m\ + 2m 2 -, -\ — . 



m% = m\_ = -(m\ + m\ + 2m 2 D + -M 2 Z cos 2/3) (31 



+ -\ Um 2 D (A u - ficot /3) 2 + (m\-m 2 - -M|cos2/3) 2 
2 V 2 



and they combine to form two complex scalar fields, 



V\ = —=(i' + + iu_) = cos 9 v L — sin^ v R , (39) 

v2 

v 2 = ~^(N+ + iN-) = sin 9 u L + cos 9 v R (40) 

a/2 



with mp 1 = mp ± , m S2 = m^ , 9 = 9 + = 9^., and 

sin 2^~ = 2m c (^-^cot/3) (41) 

Um 2 D (A u - fjt cot /3) 2 + (m| -m\- \M 2 cos 2/3) 2 

Notice that these two sneutrino states, V\^ are equivalent to the usual sfermion mass 
eigenstates within the MSSM. 

In this Dirac limit it is interesting to study the similarities in the analytical behavior of 
the neutrino/sneutrino radiative corrections and the other MSSM fermion/sfermion ra- 
diative corrections. In particular we are interested in the comparison with the top/stop 
radiative corrections. As for the phenomenological implications, this limit is not ex- 
pected to lead to relevant numerical results, since to get compatibility with the exper- 
imentally tested small neutrino masses, \m v \ ~ 0.1 — 1 eV one needs Yukawa couplings 
extremely small, Y u ~ 10~ 12 — 10~ 13 . 

(3) The MSSM limit: 

This limit is reached when one sets thd = (the value of thm is not relevant since 
once the Yukawa couplings are set to zero the predictions are absolutely independent 
of this mass scale) and one is left with a neutrino/sneutrino sector with just pure 
gauge couplings. Concretely, there are just interactions of the left-handed neutrinos 
and the 'left-handed' sneutrinos to the Z boson, exactly as in the MSSM. We are 
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interested in this limit, because we want to compare the radiative corrections from the 
neutrino/sneutrino sector within the MSSM-seesaw with those within the MSSM and 
to find the interesting regions in the new parameters of the MSSM-seesaw where the 
deviation from the MSSM result could be sizeable. 

3 Higher-order corrections to rrih 

3.1 The concept of higher order corrections in the 
Feynman-diagrammatic approach 

In the Feynman diagrammatic (FD) approach the higher-order corrected CP-even Higgs 
boson masses in the MSSM, denoted here as M^ and M# (the corresponding masses in the 
MSSM-seesaw model are denoted as M^' u and M v n ; ), are derived by finding the poles of the 
(h, H )-propagator matrix. The inverse of this matrix is given by 

(A s-i ■( p 2 -m 2 H + t HH {p 2 ) t hH {p 2 ) \ ( , 

(AHiSgs) -~\ t hH ^) p 2 -ml + Z hh ( P >)) ■ (42) 



Determining the poles of the matrix AHi ggs in (142 p is equivalent to solving the equation 

p 2 - m 2 h + ± hh (p 2 )] [p 2 -m 2 H + £hh(p 2 )] - [z hH (p 2 )] 2 = . (43) 

In perturbation theory, a (renormalized) self-energy is expanded as follows 

£(p 2 ) = £ (1 V) + £ (2 V) + --- , 

£(p 2 ) = E^V) + £ (2 V) + • • • , (44) 

in terms of the zth-order contributions E^, E^. In the following sections we concentrate on 
the one-loop corrections and drop the order index, i.e. E = E^ 1 ) in the following. 

3.2 One-loop renormalization 

In order to calculate one-loop corrections to the Higgs boson masses, the renormalized Higgs 
boson self-energies are needed. Here we follow the procedure used in [T5l[33] (and references 
therein) and review it for completeness. The parameters appearing in the Higgs potential, 
(J18)) . are renormalized as follows: 

M|->M| + 5M|, Th ^T h + 5T h , (45) 

M 2 W ^M 2 W + bM 2 Wl T H ^T H + 5T H , 

M Hi gg s -> M Hi gg s + 5M Hig g s 5 tan P -> tan/3 (1 + 5tan/3 ). 

M^j denotes the tree- level Higgs boson mass matrix given in ( !25l) . Th and Th are the 
tree-level tadpoles, i.e. the terms linear in h and H in the Higgs potential. 
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The field renormalization matrices of both Higgs multiplets can be set up symmetrically, 

h\ (l + \5Z hh \8Z hH \ (h 



- 1 *z„ lH r (46) 



Hi \ \8Z hH l + \5Z HH l \H 



For the mass counter term matrices we use the definitions 

,^2 ( 6l7l h Sm hH\ 

5M?,- = . (47) 

Hlsgs \Sml H 5m 2 H ) 

The renormalized self-energies, £(p 2 ), can now be expressed through the unrenormalized 
self-energies, S(p 2 ), the field renormalization constants and the mass counter terms. This 
reads for the CP-even part, 

± hh (p 2 ) = E hh (p 2 ) + 5Z hh (p 2 - m 2 h ) - 5m 2 h , (48a) 

± hH (p 2 ) = T, hH {p 2 ) + 5Z hH (p 2 - \{m\ + m 2 H )) - 5m 2 hH , (48b) 

±hh(p 2 ) = Z HH (p 2 ) + 5Z HH {p 2 - m 2 H ) - 5m 2 H . (48c) 

Inserting the renormalization transformation into the Higgs mass terms leads to ex- 
pressions for their counter terms which consequently depend on the other counter terms 
introduced in (145|) . 

For the CP-even part of the Higgs sectors, these counter terms are: 

5m\ = 5 Ml cos 2 (a - (3) + 5M 2 Z sin 2 (a + (3) (49a) 

+ 2M z l w c w (3Th cos(a - j3) sin 2 (a - j3) + 5T h sin(a - /3)(1 + cos 2 (a - (3))) 
+ 5tan/3 sin/3 cos/3 (M 2 A sin2(a - /3) + M§ sin2(a + /3)), 
Sm 2 hH = \{5M\ sin 2 (a - /3) - 5M 2 Z sin 2 (a + (3)) (49b) 

+ 2l7Fb^ T " sin3 (« - Z 3 ) - ^cos 3 (a - 13)) 
-5tan/3 sin/3 cos/3 (M^cos2(a - /3) + M§ cos2(a + /3)), 
5m 2 H = 5 Ml sin 2 (a - /3) + 5M Z cos 2 (a + (3) (49c) 

— 2M e s c (5Thcos(o! — /3)(1 + sin 2 (a — (3)) + 5Thsm(a — (3) cos 2 (a — (3)) 
-<5tan/3 sin/3 cos/3 (M^sin2(a - /3) + M|sin2(a + /3)) . 

For the field renormalization we choose to give each Higgs doublet one renormalization 
constant, 

Hi -)• (1 + \5Z nx )U x , U 2 ^{l + \5Z H , 2 )U 2 . (50) 

This leads to the following expressions for the various field renormalization constants in ([46]) : 



SZhh = sin 2 a SZ^ + cos 2 a 5Z-^ 2 , (51a) 

5Z hH = sin a cos a {5Z n , 2 — 5Z Hl ), (51b) 
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8Z HH = cos 2 a8Z Hl + sm 2 a8Z H2 . (51c) 

The counter term for tan (3 can be expressed in terms of the vacuum expectation values as 

5 tan/3 = \ (5Z H2 - 8Z Hl ) + ** - ** , (52) 

2 v 2 V\ 

where the 8vi are the renormalization constants of the vf. 

ui ->(1 + SZ m ) ( Vl + S Vl ) , v 2 -)■ (1 + <5Z % ) ( U2 + <fo 2 ) . (53) 

It can be shown that the divergent parts of 5vi/vi and 8v 2 /v 2 are equal [15] . Consequently, 
one can set 5v 2 /v 2 — 8v\/v\ to zero. 

The renormalization conditions are fixed by an appropriate renormalization scheme. For 
the mass counter terms on-shell conditions are used, leading to: 

8M 2 z = ReZ zz (M 2 z ), 8M 2 w = ReT, ww {M 2 w ), 8M 2 A = ReZ AA (M 2 A ). (54) 

Here T>zz,ww denotes the transverse part of the self-energies. Since the tadpole coefficients 
are chosen to vanish in all orders, their counter terms follow from Ti^m + STi^m = 0: 

8T h = -T h , 8T H = -T H . (55) 

For the remaining renormalization constants for 8 tan (3, 8Z-^ 1 and 8Z-u 2 various renormal- 
ization schemes are possible [3TH33] . 

On-shell renormalization 

One possible choice is an on-shell (OS) renormalization. The renormalization conditions for 
the renormalized Higgs-boson self-energies are 

(56) 

(57) 



(58) 
(59) 











K h «) = o , 










sWK) = o . 


This yields 








8Z™ = -ReZ> hh (m 2 h ), 
8Z° s H = -ReX' HH (m 2 H ), 


equivalently 


to 










6Z% = 


1 
cos 2a 


(sin 2 a Re'E' hh (m 2 l ) — cos 2 a ReE' HH (?7 




8Z™ = 


1 
cos 2a 


(" 


- cos 2 a Re 'E' hh (m 2 l ) + sin 2 a Re T! HH 


For 8 tan /3° £ 


1 a convenient choice is 
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)) , (60) 



m 



li 



)) • (61) 



= j^^ (^ Kk(m 2 h ) - Re Tf HH {m%)) . (62) 

It should be kept in mind that this scheme can lead to large corrections to m>h in the 
MSSM [3TJIM] . hence worsening the convergence of the perturbative expansion. Furthermore, 
it is known to provide gauge dependent corrections at the one-loop level 



13 



DR renormalization 

A convenient choice which avoids the previously commented large corrections to rrih in the 
MSSM and is (linear) gauge independent at the one-loop level is a DR renormalization of 
<5tan/3, 5Z Hl and SZ n2 , 

5Z^=-[ReE' HHla=0 f\ (63a) 

5Zgf=-[ReS;, |Q=0 ] div , (63b) 

5t,n^= 1 -(5Z™-5Z™) . (63c) 

The [ ] lv terms are the ones proportional to A = 2/e — 7E + log(47r), when using dimensional 
regularization/reduction in d = 4— e dimensions; 7e is the Euler constant. The corresponding 
renormalization scale, /ipr^, has to be fixed to a certain mass scale that will be discussed 
below. 



Modified DR renormalization (mDR) 

The fijyK dependence introduced in the DR scheme can lead in the present context to large 
logarithmic corrections oc log(rra| / //i^-) for large values of the Majorana mass tum (as will 
be discussed below). These large corrections could again worsen the convergence of the 
perturbative expansion. One possible way out is to replace [ ] lv by [ ] m lv , where the latter 
means to select not only the terms oc A as in (1631) . but the terms oc A m = A — log (m 2 M /fi^). 
This prescription for the counterterms defines the modified DR renormalization scheme, 
which will be named in this work in short as mDR, 

^^ = -[ReSU,^o] mdlV , (64a) 

5Z^ = -[ReJ:' hh ^ ] mdi \ (64b) 

5tan/3 mm = - (6Z%™ - 8Z^\ . (64c) 

As will be shown below, effectively this corresponds to the particular choice of hym = m M- in 
this way the potentially large logarithms vanish, what makes it a convenient choice. Usually 
this choice is referred to in the literature as 'decoupling the large mass scale by hand' (see 
e.g. [35j[36] and references therein). 

It should be kept in mind that in the mDR scheme the parameter tan (3 = tan / g mDR has 
a different meaning than the "conventional" parameter tan/3 = tan/3 DR . However, we have 
checked that this shift is numerically insignificant. 

4 Results 

In this section we first present the results of the one-loop corrections from neutrino/sneutrino 
contributions to the neutral Higgs boson renormalized self-energies within the MSSM-seesaw 
and then we discuss the derived results for the Higgs mass corrections. 
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4.1 One-loop calculation of the renormalized self-energies 

The full one-loop neutrino/sneutrino corrections to the self-energies, £^\ £## and £^, 
entering (1431) have been evaluated with the help of FeynArts I37) 1 ! and FormCalc [38]. For 
shortness, in this and the next subsection these self-energies will be named simply as E^, 
T'HH, and T,hH, respectively. The new Feynman rules for the neutrino/sneutrino sector, 
derived in this work and collected in the Appendix A, have been inserted into a new model 
filezl As regularization scheme we have used dimensional reduction [39], thus preserving 

SUSY gnum. 

The generic one-loop Feynman-diagrams contributing to the renormalized self-energies 
are depicted in Fig. [TJ They include the two-point and one-point diagrams in the Higgs 
self-energies, tadpole diagrams, and the two-point and one-point diagrams in the Z boson 
self-energy Here the notation is: <j) refers generically to all neutral Higgs bosons, h, H, A; F 
refers to all neutrinos n.j (i = 1, 2); S refers to all sneutrinos n« (i = 1, ..4), and Z refers to 
the Z boson. 




F 




S 



F 



: s 



■--.- .' N > s 



1 , ' I 



I fh \ 







I 



\ r ^- . 



\ 



r ■ 



s 








Figure 1: Generic one-loop Feynman-diagrams contributing to the neutral Higgs bosons 
renormalized self-energies (see text) 

The analytical results for the unrenormalized self-energies and tadpoles are collected in 
the Appendix B. The final analytical results for the renormalized self-energies are easily 
obtained by inserting these results into ( 1481) . 



lr rhc program and the user's guide are available via www.feynarts.de. 
2 This model file is available upon request. 
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We have checked that all the divergences involved in the computation cancel and the 
renormalized self-energies, E^(p 2 ), E##(p 2 ) and E/^p 2 ) in the three schemes OS, DR, 
and mDR are all finite, as expected. We have also checked that the renormalized self- 
energies in the OS scheme, are independent of the regularization scale /igR, as they must be. 
The renormalized self-energies in the DR are /ijjg; dependent whereas the ones in the mDR 
scheme are /x^ independent by construction. Analytically they are related by E mDR (p 2 ) = 

fpR( n 2\\ 

4.2 Analysis of the renormalized self-energies 

In the following we discuss the numerical results for the renormalized self-energies. They are 
collected in Figs. l21through[T(Il First we compare the predictions of the one-loop renormalized 
self-energies in the three schemes for the full interval 10 3 GeV < tum < 10 15 GeV, and 
next we analyze these exact results at large tjim with the help of the simple analytical 
formulas that are obtained in the seesaw limit. Then we choose the mDR scheme and 
show the exact numerical results of the renormalized self-energies as functions of all the 
neutrino/sneutrino parameters involved. Finally we conclude on the subset of most relevant 
parameters (specifically, m M , m^, B u and m„) which will be the selected ones to study the 
corrections to M^ in the next subsection. For the final estimate of these corrections, and 
to localize the regions of the parameter space where they can reach sizeable values, we will 
vary these relevant parameters within some selected plausible intervals. For the parameters 
which do not exhibit a relevant numerical effect on Mh (specifically, tan/3, Ma, (J>, mi and 
A v ) we choose representative values. For completeness, we will also comment shortly at the 
end of this subsection on the Dirac case. 

In order to compare systematically our predictions of the neutrino/sneutrino sector in the 
MSSM-seesaw with those in the MSSM, we have split the full one-loop neutrino/sneutrino 
result into two parts: 

£(P 2 ) | full = S(p 2 ) | gaugc + t(p 2 ) | Y ukawa , (65) 

where E(p 2 ) | gaU ge means the contributions from pure gauge interactions and they are obtained 
by switching off the Yukawa interactions, i.e. by setting Y v = (or equivalently mjj = 0). 
The remaining part is named here E(p 2 )| Yu kawa an d refers to the contributions that are only 
present if Y v ^ 0. In other words, this separation splits the full result into the common 
part with the MSSM, given by E(p 2 )| gaugc , and the new contributions due to the presence of 
Majorana neutrinos with non vanishing Yukawa interactions, given by S(p 2 )| Yukawa- Thus, 
by comparing the size of these two parts, within the allowed parameter space region, we will 
localize the areas where E(p 2 )| Yukawa ~^> E(p 2 )| gaugc , which will therefore indicate a significant 
departure from the MSSM result. 

Dependence on tum 

We show in Fig. [2] the predictions for E^Qo 2 ) as a function of vtiM i n the three schemes: 
DR (upper left plot), OS (upper right plot), and mDR (lower left plot). In these plots we 
have considered an extremely wide range for the tum values, from 10 3 GeV up to 10 15 GeV, 
and fixed the physical light neutrino mass to \m v \ = 0.5 eV. Consequently, nip is derived 
from tum and m v by using (ITT1) and (fT2|) . The other parameters are fixed as indicated in 
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the figure. In this and in the following figures we have fixed p 2 in the self-energies to a 
particular value, corresponding to an approximation of the higher-order corrected value of 
Mh for the input MSSM parameters set in each figure, see below. The numerical values used 
here and in the following for the SUSY parameters are representative values (as will also be 
shown below). Therefore, despite choosing only a few values for the parameters, the results 
obtained can be considered as more general. 
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Figure 2: Renormalized Higgs boson self-energies as a function of ttlm and comparison 
between the three considered schemes. Upper left panel: S^(p 2 ). Upper right panel: 
t^(p 2 ). Lower left panel: t^(p 2 ). Lower right panel: t^(p 2 ), t^(p 2 ) and tf^(p 2 ). 
All self-energies are evaluated at p 2 = (116 GeV) 2 . 



In the three mentioned plots in Fig. [2] one can see that the numerical value of the full 
result is nearly constant with tum in the three schemes from vim = 10 3 GeV up to thm ~ 
10 12 GeV. Furthermore, this constant value is approximately the same in the three schemes 
(the differences are below ~ 1CT 2 GeV 2 ), and is totally dominated by the 'pure gauge 
contributions'. Thus, for 10 3 GeV < ttim ~ 10 12 GeV the result in the MSSM-seesaw nearly 
coincides with the result in the MSSM, irrespectively of the scheme. For the choice of input 
parameters in this plot, we get f^hh\faii — S^aI™.™ — —23.67 GeV 2 . 



J hh\ gauge 
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For larger values of vtiu in the range 10 12 GeV < tjim < 10 15 GeV, there are, however, 
remarkable differences between the three considered schemes, and the main differences come 
clearly from the 'Yukawa contributions'. Whereas £^|f u ii is apparently constant with ttim, 
also for iriM > 10 12 GeV, |£° R |f u n| and |S^ DR |f ull | grow noticeably with m M at these large 
vfiu values. The numerical value of S^|f u n is negative for m M > 10 12 GeV and gets large 
values in this range, where they are totally dominated by the 'Yukawa contributions'. For 
instance, for m M = 10 13 GeV, we get £° R |f u n ~ £° R | Yukawa — —250 GeV 2 , and for m M = 
10 14 GeV, we get iffifai ^ Sj^lvukawa ^ -3000 GeV 2 . In the mDR scheme, the result is 
negative up to 5 x 10 13 GeV and then becomes positive and large for ttim > 5 x 10 13 GeV. 
Notice that, the absolute value in the mDR scheme at large ttim is always smaller than in 
the DR scheme, due to the commented cancellation of the large logarithms \og(rriM/ ^Ym) 
corresponding to the choice /i^ = rriM- Notice also that, in spite of this cancellation, the 
size of the corrections in mDR, are still large for large enough tum values. For instance, for 
m M = 10 15 GeV, we get dominance of the 'Yukawa contributions' S^ DR |f u ii ~ £™/P R |Yukawa — 
500 GeV 2 . In contrast, for tum = 10 14 GeV, the Yukawa contributions' and the 'pure gauge 
contributions', compete since S^ DR | Yu kawa — 60 GeV 2 and S™ / P R | gaugc ~ —24 GeV 2 leading 
to E£pkn ~ 36 GeV 2 . _ 

In the lower right plot of Fig. [2] we compare E^ J - ,R |f u n to the other two renormalized 
self-energies, £^§- R |f u n and S^ R | full . One can observe that the three self-energies behave 
qualitatively very similarly with m,M, being approximately constant for rriM < 10 12 GeV 
and growing (in modulus) with itlm for 10 12 GeV < vtiM < 10 15 GeV. For the choice of 
parameters in this plot, |E^P R |f u n| is larger than the others in the full explored niu range. 
This will be relevant for the forthcoming estimate of the one-loop radiative corrections to M^. 

The previously commented growing behavior of the renormalized self-energies with tum 
is a consequence of the corresponding growing behavior of the neutrino Yukawa interactions 
with rriM, see fTTTT) and ( TT2l) . This is a well known feature of the seesaw models that, in order 
to get the light neutrino masses m v in agreement with data, one must impose for each input 
vtiM value the proper Y v (and therefore mo) to precisely match the experimentally inspired 
input m u . Y v is therefore not an input but an output in this approach, and according to ffTTT) 
and ( [121) Y u grows with tum as Y v oc y/rriM- The behavior of the renormalized self-energies 
with rriM is, consequently, the result of the two competing facts, the increase of Y v with tum 
and the decreasing with tum from the neutrino and sneutrino propagators in the loops. 

Dependence on tum in the seesaw limit 

In order to illustrate more clearly the behavior with thm, we have analyzed in more detail 
the renormalized self-energies in the seesaw limit, as defined in section [2j As the increase 
with jum starts at very large ttim > 10 12 GeV values (i.e. much larger than the other scales, 
m M 3> m D , M z , Ma, fi, 1711,171^, B U ,A U ), one expects that this limit should approximate 
pretty well the full result and show its same main features. 

For the computation of the renormalized self-energies in this seesaw limit, we have 
performed a systematic expansion of the exact result in powers of the seesaw parameter 
£ = thd/itim- In order to reduce the number of parameters, and for a clearer interpretation 
of the results, we have set in this expansion, A v = fi = B v = (which is justified, see below) 
and we have assumed universal soft SUSY breaking masses, i.e., m^ = m^ = ttisusy- 



The analytical expressions for these expanded renormalized self-energies are of the generic 
form: 
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is the first term in the expansion, i.e. O(£ ), (S(p 2 )) is the next 
S(p 2 ) ) is the term of 0(£ 4 ), etc. It should be noticed that there are 
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no terms with odd powers of £. The first term in this expansion is precisely the pure gauge 
contribution, (S(p 2 ) J = S(p 2 )| gauge . Therefore, it approximates the result in the MSSM 

and the rest approximates the Yukawa part, 
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In order to get simple formulas, we have expanded in addition each term in the series in (1661) 
in powers of the other small dimensionless parameters, namely, Mz/rriM, Ma/trm, p/ttim 
and m SVS Y/m M . 

The result of the previous seesaw expansion (we just show the leading terms; terms 
suppressed by factors \jvr? M respect to these leading ones are not relevant and, therefore, 
are not included) for each of the three considered renormalization schemes is as follows. 
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Figure 3: Comparison between the predictions from the seesaw expansion and the exact 
results for the Yukawa part. Left panel: mDR scheme. Right panel: OS scheme. In both 
panels, p 2 = (116 GeV) 2 . 



From these formulas the qualitatively different behavior of the renormalized Higgs-boson 
self-energies on the Majorana mass scale itim can be understood. The main difference be- 
tween the OS scheme and the DR/mDR schemes appears in the Yukawa part, especially in 
the term of 0(m 2 D ). At the various orders the comparison of the three schemes is given as 
follows. 

At the leading order in the seesaw expansion, O(m D ) in (|68l) . the results in the DR and 
mDR schemes coincide. This is indeed a consequence of the fact that, at this order, E^(p 2 ) 
turns out to be /a^ independent. The result in the OS scheme differs from these later by a 
term of order (yf 2 M|M| w /77ig USY , where M^ w refers generically to the involved masses of the 
order of the electroweak scale, i.e., M|, p 2 , M§, m 2 tree . Furthermore, this difference turns 
out to be numerically extremely small. This explains why, for low values of the Majorana 
scale, where the O(m D ) term of the expansion dominates, the predictions from the three 
schemes are nearly indistinguishable. 

At the next order in the seesaw expansion, 0{m 2 D ) in (155J1 . the OS result differs substan- 
tially from the DR and mDR schemes. First, the OS result is extremely suppressed with 
respect to the DR and mDR results at large mM- This is due to the fact that the leading 
contribution, i.e. of the order of g 2 m 2 D M'^ N jM\ 1 vanishes in the OS whereas it is present in 
the other schemes. As can be seen in (169 p . the first non vanishing contribution contains an 
extra factor ~ m susY/ m M which can be extremely small for tum ^ m susY- This remark- 
able difference of the OS result has its origin in the different values of the SZhh and S tan /3 
counterterms. More specifically, by computing their finite parts in the OS scheme and in the 



21 



seesaw limit, we get 



5 os Z 



hh | finite 



5 Ub tan/3| finite 



g 2 m 2 D cos 2 a 
' 32c 2 w M 2 z -k 2 sin 2 (3 

2 2 
9 m D 

' Uc 2 M 2 tt 2 sin 2 



log 



■//? 



M 



/'■ 



DR 



log 



??? 



M 



[I- 



DR 



+ o 



+ C 



^EW' m SUSY 



■/?? 



iW 



Ml w , 



??? 



SUSY 



m 



M 



(71) 
(72) 



These finite contributions lead to the cancellation of the above commented leading contri- 
butions. 

In the DR scheme, we get an explicit logarithmic dependence on rriM, concretely as 
— \og{rn 2 M / jji 2 ^) . By construction this term is absent in the mDR result. Therefore, the 
main difference between these two schemes DR and mDR is this logarithmic contribution 
that can be sizeable for very large thm 3> /^dr- 

The results at the next to next order in the seesaw expansion, 0(raf)) in (ITU]) , show that 
they all go (leaving apart the logarithms) as g 2 m^(M! w ,mg USY )/(M§m 2 . f ). Therefore the 
0(m%) terms are extremely suppressed in the three schemes, and consequently they are not 
relevant in the large tum regime. 

All the above commented analytical features of the seesaw expansion have also been 
checked numerically, as it is illustrated in Fig. |3j In this figure we show separately the 
0(m 2 D ) and 0(m%) contributions and the exact Yukawa prediction in both the mDR (left 
plot) and OS scheme (right plot)! 3 ] One clearly observes the dominance of the 0(m 2 D ) over 
the 0(m%) in the mDR scheme by many orders of magnitude in the full explored m M range. 
One also sees that the 0{m 2 D ) result approximates extremely well the exact Yukawa result for 
rriM ^ 10 4 GeV. In contrast, in the OS scheme, the 0{m 2 D ) term dominates just up to about 
rriM = 10 10 GeV, but then for larger values the 0(m%) dominates. In this plot it is also 
manifested that the exact Yukawa result in the OS is well approximated by the (Dim 2 -)) term 
in the interval 10 3 GeV < m M < 10 11 GeV and by the 0(m%) term for m M > 10 12 GeV. At 
this large values, however, the size of the correction is extremely small (below 1CT 17 GeV 2 ), 
hence, irrelevant. It is also clear from this plot that the numerical results for the 0(712%) 
contributions are similar in the three schemes. 

From the definition of the three renormalization schemes, see Sect. 13.21 and our analytical 
and numerical analysis in this section we conclude that the mDR scheme is best suited 
for higher-order calculations in MSSM-seesaw model. The other two schemes can lead to 
unphysically large corrections at the one-loop level. We will focus in the following on this 
scheme, and the numerical evaluation of M^ 



see Sect. 14. 3[ will be performed solely in this 
"preferred" scheme. 

Finally, in this context, we discuss the decoupling or non-decoupling behavior of the 
neutrino/sneutrino one-loop radiative corrections with the Majorana scale. According to 
Figs. |2] and |3l the Yukawa part of the renormalized self-energy in the mDR scheme grows 
with Vfiw However, this does not constitute by itself a proof of non-decoupling of tum in 



3 It should be kept in mind that due to the different renormalization of tan j3 the meaning of this input 
parameter is different in OS and in the mDR scheme. In order to perform a real numerical comparison a 
transition from tan/3 = tan/3 mDR — > tan/3 os would have to be performed. However, here we are interested 
in the qualitative behavior and we do not consider this shift. 
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Figure 4: Left panel: Decoupling/Non-decoupling behavior of the one-loop neu- 
trino/sneutrino corrections to the renormalized lightest Higgs boson self-energy at large 
rriM in the mDR scheme. Right panel: Dependence of the neutrino Yukawa coupling (and 
ttid) with ttim- 



the radiative corrections to S™^ 11 for asymptotically large mu- To analyze this question, we 
have to investigate separately the behaviors of £™/P R and mo with tjim, since in the way the 
seesaw mechanism is implemented here, as we have mentioned before, mjj (or equivalently 
Y v ) is not an input but an output and it grows proportional to ^m M . To analyze these two 
behaviors separately we show in the left plot of Fig. H]the ratio (Sj^P^Yukawa/^z) versus tjim 
(and mo), and in the right plot we show the predictions of the Yukawa coupling (and mo) as 
a function of m M . The latter one exhibits the (trivial) result of Y v oc ^mu as expected. In 



the left plot a constant behavior of the ratio (Sj^P^Yukawa/^l) is clearly manifested, which 
means that the growing of (S™° R ) Yu kawa with m M is exclusively due to the growing of Y v 
(or m £>) with mu- However, still this ratio turns out to be no n- vanishing for asymptotically 
large mM, arid constant with mo, as can be seen in Fig. HJ Therefore, a non- decoupling 
constant behavior must be concluded in the Majorana case from all this discussion. This 
constant, on the other hand, is very well approximated by the coefficient multiplying the 
factor m D in the S^f 1 ^ 2 )^ result of fl69|) . 

In order to understand this issue better, we compare this analytical result, showing a 
constant behaviour of the renormalized Higgs boson self-energy in the tjim -^ oo limit when 
Y v is kept fixed, with the corresponding result in the Dirac case. For simplification in this 
analytical comparison we focus just on the 0(p 2 m 2 D ) terms and use the electroweak basis for 
neutrinos and sneutrinogj. The results at 0{p 2 m D ) for the renormalized self-energies in the 
DR scheme for the Majorana and Dirac cases are: 



-\Majorana,DR 
J hh 



(P 2 ) 



q 2 m 2 n p 2 cos 2 a , 1 , m 



2 
M 



32vr 2 M 2 sin 2 /T2 



DR 



4 The computation in this case reduces to just the evaluation of one type of loop diagrams, the sunset 
diagrams, 2nd and 5th in Fig. [1] 
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g 2 m 2 D p 2 cos 2 a 

64vr 2 M2,sin 2 /3 l ' 

2 2 2 2 2 

V^Dirac,DR/ 2\ 9 ™dP cos a / n i P \ /-.s 

E » (P) = 32SP^ ,2 - |0e 4» (74) 

where the first and second lines in £^ a J° rana > (p 2 ) are the contributions from neutrinos 
and sneutrinos respectively. It should be noticed that the 0(p 2 m 2 D ) sneutrino contributions 
come exclusively from the new couplings g' h $ Li> = - ^^j^i which are not present in 
the Dirac case. It should also be noticed that this result in the Majorana case translates 
into our 0(p 2 m 2 D ) term in (169a)) . The comparison of the two formulas shows that the result 
of the Majorana case for low momenta, p 2 <C rn 2 M , does not coincide with the result of the 
Dirac case. 

From the right plot in Fig. H] we can also conclude on the range of tum values where 
the neutrino Yukawa couplings get too large and potentialy non-perturbative. The con- 
crete crossing line to set the perturbativity region is not uniquely defined, but it should 
be considered around Y v ~ Oil). For instance, by setting the crossing at Y 2 /(An) = 1.5 
(Y u = 4.34) we get perturbativity for m M < 10 15 GeV, and by setting it at Y v = 1.5 it is 
got for m,M < 10 14 GeV. In the following of this subsection we set mu = 10 14 GeV as our 
reference value. 

Dependence on tan/3, Ma, fJ>, rni, m^, A„, m v , B v and p 

The behavior of the renormalized self-energy in the mDR scheme with the other parameters 
entering in this computation are shown in Figs. [5] - [HU In all these plots we have included 
separately the gauge, Yukawa and total results for comparison. 

First, the behavior with tan/3 is analyzed in the left plot of Fig. [5j It exhibits basi- 
cally the expected features that can be inferred from the loop corrections of an up-type 
fermion/sfermion. The neutrino/sneutrino one-loop radiative corrections reach their maxi- 
mum value at the lowest considered value of tan/3, tan/3 = 2 in this plot. For tan/3 > 5 the 
dependence is nearly flat. There are no relevant differences between the behaviors with tan j3 
of the Yukawa and the gauge parts. From now on, we will set tan /3 = 5 as our reference 
value. 

The behavior with Ma is displayed in the right panel of FigJSJ Again we see no relevant 
differences with respect to the well known behavior in the MSSM. For Ma larger that 150 
GeV the total contribution from the neutrino/sneutrino sector to the renormalized self- 
energy is nearly flat with Ma- In the following we will take Ma = 200 GeV as our reference 
value. 

The dependence with the soft SUSY breaking mass of the 'left handed' SU(2) doublet, 
mi, is shown in Fig. We see that the gauge contribution is negative and increases in 
modulus with increasing m^, whereas the Yukawa contribution is positive and nearly insen- 
sitive to changes of m^ in the investigated interval, 10 2 GeV < m^ < 10 4 GeV. The total 
neutrino/sneutrino corrections, at these selected values of the model parameters, are positive 
and decreasing with m^ for 10 2 GeV < m^ < 2 x 10 3 GeV and then become negative and 
increasing in modulus with mi for 2 x 10 3 GeV < m^ < 10 4 GeV. 

The behavior with the soft SUSY breaking parameter of the 'right handed' sector m^ 
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Figure 6: E™^ R (p 2 ) as a function of m^; we have set p 2 = (105 GeV)' 



is shown in Fig. [71 In the left plot a mass scale similar to the other soft SUSY-breaking 
parameters is investigated, whereas in the right plot values of m^ closer to tum are explored. 
It should be reminded that these values are not constrained by data. An interesting feature 
can be observed at large values of m^. The contributions to the renormalized self-energy 
stay flat up to about m^ ~ 10 13 GeV. Above this mass scale the Yukawa part grows rapidly, 
reaching very large values at m^ ~ 10 14 GeV of around S™^ 11 ~ 7000 GeV 2 . 

The behavior with the new soft SUSY-breaking trilinear coupling A v is shown in the 
left plot of Fig. [8] The full result, the gauge, and Yukawa parts are nearly independent on 
this parameter in the studied interval, —1000 GeV < A v < 1000 GeV. Although not shown 
explicitly, we have also studied the behavior with /i and got the same 'flat' behavior for 
— 1000 GeV < /i < 1000 GeV. This justifies our choice A u = fi = in our seesaw expansion 
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The behavior with the lightest neutrino mass, 



m. 



is demonstrated in the right plot of 



Fig. [8] One can see that the Yukawa part is quite sensitive to this mass that we have varied 
in a plausible and compatible with data range. The growing of the result with \m u \, for fixed 
mui is the consequence of the growing of Y v (or mn) with \m v \ since in this model they are 
correlated, as shown in ( ITTl) and ( TT2l) . 

The behavior with B v is analyzed in Fig. |9j We have found a flat result with this new 
soft parameter for most of the explored range, except at very large values, B v > 10 12 GeV, 
as shown in the right plot. For these large values the Yukawa part grows noticeably with B u 
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Figure 9: E™° R (p 2 ) as a function of B v . Left panel: low B v values, 10 2 GeV < B v < 
10 4 GeV. Right panel: high B v values, 10 12 GeV < B u < 10 13 GeV. In both panels we have 
set p 2 = (105 GeV) 2 . 



and dominates largely the total result, leading to large radiative corrections. For instance, 
for the parameters chosen in this figure and B v = 10 13 GeV, we found £™/P R ~ 2400 GeV 2 . 
The question whether such large values of B v are realistic depends on the particular models 
and universality conditions. However, such an analysis is beyond the scope of our paper. 
On the other hand, if we apply the bounds that are imposed in order to avoid destabilizing 
the electroweak symmetry breaking [13], leading to B U Y 2 / (Sn 2 ) < mgusy/tan/S, one gets 
an upper limit on B u . For Y u ~ 1, m SUSY ~ 1000 GeV and tan/3 ~ 5 one finds B u < 
1.6 x 10 4 GeV. For this range the renormalized Higgs-boson self-energy is nearly independent 
of B v . From now on, we will choose B v = 500 GeV as our reference value. 

Finally, we show in Fig. [10] the behavior with p 2 , the square of the external momentum 
of the Higgs boson self-energies, which is a relevant issue for the discussion of the radiative 
corrections to the Higgs-boson masses (see the next subsection). The three renormalized 
self-energies, E^, ^hh and XW> are clearly dependent on p 2 , but the most sensitive one is 
S/ift. It is clear from this figure that setting p 2 = in the renormalized self-energies does 
not provide a good approximation for the estimate of the radiative corrections to the Higgs 
boson mass from the neutrino/sneutrino sector in the present case of Majorana neutrinos. 
One can also see that mainly the Yukawa part is responsible for this sensitivity to p 2 . Setting 
the proper p 2 in order to estimate realistically the Higgs mass corrections will be discussed 
in the next subsection. 



The Dirac case 

Finally, we perform a comparison between the case of massive Majorana neutrinos (as an- 
alyzed so far) and the case of Dirac neutrinos. In order to analyze the Dirac case, we have 
computed the one-loop neutrino/sneutrino contributions to the renormalized lightest Higgs 
boson self-energy for thm — 0. The analytical results for this Dirac case are collected in 
Appendix C. We have chosen here the DR scheme, since due to the absence of niu no large 
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logarithmic corrections are expected, and a comparison to existing calculations can readily 
be performed. First, we have checked the fmiteness of the result. Second, we have also 
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Figure 11: One- loop corrections to the Yukawa part of the lightest Higgs boson renormalized 
self-energy from the neutrino/sneutrino sector in the case of massive Dirac neutrinos 



checked that the obtained formulas agree with the well known result of the one-loop radia- 
tive corrections from other massive fermion/sfermion sectors of the MSSM, with the obvious 
corresponding changes of fermion/sfermion parameters and quantum numbers. In particu- 
lar, it can be seen that the formulas in Appendix C coincide with the one-loop corrections 
from the MSSM top/stop sector by replacing, correspondingly, the neutrino SU(2) x U(l) 
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As for the numerical estimate, we present in FigfTJJthe result of the Yukawa contributions 
from the one-loop neutrino/sneutrino radiative corrections to the renormalized self-energy, 
(S^ kawa )Dirac 5 as a function of the physical neutrino mass, \m v \ = rri£>. The regularization 
scale has been fixed here to /^g^ = 100 GeV and the external momentum to p = 116 GeV. 
As in the Majorana case, we consider an interval for the neutrino mass inspired by exper- 
imental data, 0.01 eV < \m„\ < 1 eV. In this plot we see clearly that, as expected, these 
Yukawa contributions are extremely small (below 10~ 20 GeV 2 ) and are fully dominated by 
the gauge part which we have also estimated, for the chosen parameters in this plot, leading 
to (S|^ ugc )Dirac = —18.5 GeV 2 . Notice that this gauge part is similar in both Majorana and 
Dirac cases, as can be seen in the right plot of FigJHl In summary, the radiative corrections 
from the massive neutrinos/sneutrinos in the Dirac case are phenomenologically irrelevant 
and therefore this case is totally indistinguishable from the MSSM with massless neutrinos. 

4.3 Estimate of the one-loop corrections from neutrino/sneutrino 
sector to Mh within the MSSM-seesaw 

We recall that the anticipated LHC precision of the mass of an SM-like Higgs boson is 
~ 200 MeV, and that at the ILC an accuracy on the mass could reach the 50 MeV level. 
These experimental precisions set the goal for the theoretical accuracies. 

As outlined in Sect. 13.11 the higher-order corrected light MSSM Higgs-boson mass is 
obtained as a pole from (|43p . i.e. where p 2 = M%. A realistic evaluation requires to take into 
account all known higher-order corrections to the renormalized Higgs-boson self-energies [32]- 
In order to simplify our analysis, but to maintain the high accuracy we follow a slightly 
different strategy. For a given set of SUSY parameters we first calculate Mh and Mh in the 
MSSM with the help of FeynHiggs [THl E31I33J 133]- In this way all relevant known higher- 
order corrections are included, but no vjv contributions are taken into account yet. This 
corresponds to a 'diagonalization' of the CP-even Higgs sector in the MSSM without heavy 
Majorana (s)neutrinos. In a second step we search for the poles of 
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where, ^hhHHhH denote the full corrections to the renormalized Higgs-boson self-energies 
from the vjv sector, obtained in the mDR scheme as described in the present work. The 
pole, the light Higgs mass including the vjv corrections (i.e. in the MSSM-seesaw model), 
is denoted by M v h . This 're-diagonalization' now effectively takes into account the full 
result of the MSSM-seesaw. The momentum in the self-energies is fixed to the value Mh 
as obtained with FeynHiggs, since it is expected that the new contributions only give a 
relatively small correction to this Mh- In a more elaborate analysis the renormalized self- 
energies should be evaluated with free p 2 . However, we expect only a very minor effect 
from fixing the external momentum to this value. In the near future the results of the new 
neutrino/sneutrino corrections will be implemented into the code FeynHiggs. 

The numerical results for Am™ DR := M v h — Mh are summarized in Figs. IT21 through [T51 
We have chosen here to explore the Higgs mass predictions as a function of just the most 
relevant model parameters which, according to our previous exhaustive analysis of the renor- 
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Figure 12: One-loop corrections to the lightest Higgs boson mass from the neutrino/sneutrino 
sector as a function of the heavy Majorana mass for various choices of the soft mass m^. 
Left panel: m R < 10 13 GeV. Right panel: 10 13 GeV < m k < 10 14 GeV. 



malized Higgs-boson self-energies, are going to provide the most interesting/sizeable correc- 
tions. These are: the Majorana mass vnu (or, equivalently, the heaviest physical Majorana 
neutrino mass ttin), the soft SUSY breaking parameters m^ and B v and the lightest physical 
Majorana neutrino mass m„. As for the numerical values of these relevant parameters, we 
focus here in the following intervals: 10 13 GeV < m^ < 10 15 GeV, 0.1 eV < \m u \ < 1 eV, 
10 3 GeV < m^ < uim and 10 3 GeV < B v < 4 x 10 12 GeV. For the remaining model param- 
eters, tan/3, M A , fj,, m^ and A u , we choose here the same reference values as in the previous 
subsection. The corresponding predictions for other choices of the parameters can be easily 
inferred from our previous results of the renormalized self-energies. 

In Fig. [12] we show the predictions for Am™ DR as a function of the Majorana mass tum, 
for several input m^ values. As a general feature, the Higgs mass corrections for the reference 
parameter values in the left plot are positive and below 0.1 GeV if tum ~ 5 x 10 13 GeV and 
mp, < 10 12 GeV. For larger Majorana mass values, the corrections get negative and grow 



up to a few GeV. For instance, Am" 



-2.15 GeV for m M 



- m M = 10 15 GeV. The results in 
ma > 10 13 GeV the Higgs mass 
corrections are negative and can be sizeable, a few tens of GeV, reaching their maximum 



the right plot show that for larger values of the soft mass, 



values at m R ~ 

Am mDR 



iriM- For instance, for m^ = vhm = 10 14 GeV we get a very large correction, 



—50 GeV. This last large negative value is in agreement with the prediction 
in Ref. [12] for the same corresponding input values of the parameters in their split SUSY 
scenario. It should be noticed that, in the case of such large corrections our approximation 
of ( |75|) is not accurate enough to obtain a precise result for M v h . However, our method still 
yields an indication of the size of the corrections from the v jv sector to M^. 

The behavior of the Higgs mass corrections as a function of the B v parameter is displayed 
in the left plot of Fig. [13j Again, Am™ DR gets negative and large for large B u , reaching the 
maximum size at B v ~ vnu ■ For instance, for the input model parameters in this plot, and 
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Figure 13: Left panel: One-loop corrections to the lightest Higgs boson mass from the 
neutrino/sneutrino sector as a function of the heavy Majorana mass, ttim, for various choices 
of the soft 5-parameter, 10 3 GeV < B v < 4 x 10 12 GeV. Right panel: Dependence of the 
Higgs mass corrections with the lightest neutrino mass, \m v \. 



The dependence of the mass corrections with the light Majorana neutrino mass is illus- 
trated in the right panel of Fig. [13j The size of the corrections grow with \m u \, as expected, 
and can be either positive in the low region, close to \m v \ ~ 0.1 eV, or negative in the high 
region, close to \m„\ ~ 1 eV. 

These same interesting features of the Higgs mass corrections in terms of the two relevant 
physical Majorana neutrino masses, m^ and m v , are summarized in the contour-plot in 
Fig. dU Here we have fixed all the soft parameters, including m^, to be at 1 TeV. The 
contour-lines for fixed Am™ DR range from positive values around 0.1 GeV in the left lower 
corner of the plot, corresponding to neutrino mass values of \m u \ =0.1 — 0.3 eV and m^ = 
3 x 10 13 GeV, up to negative values around —5 GeV in the right upper corner of the plot, 
corresponding to, for instance, \m v \ = 1 eV and m^ = 10 15 GeV. It should be noticed 
that the contour-line with fixed Am™ DR = 0.09 (drawn with a wider black line in this plot) 
coincides with the prediction for the case where just the gauge part in the self-energies have 
been included. This means that 'the distance' of any other contour-line respect to this line 
represents the difference in the radiative corrections respect to the MSSM prediction. 

We plot in Fig. [TSl the contour-lines for fixed Am™ DR in the less conservative case where 
rrift is close to m M . These are displayed as a function of \m„\ and the ratio m^/m M . m M 
is fixed here to the reference value, mu = 10 14 GeV. For the interval studied here, we see 
again that the radiative corrections can be negative and as large as tens of GeV in the upper 
right corner of the plot. For instance, Am™ DR = —30 GeV for m M = 10 14 GeV, \m v \ = 0.6 
eV and m^/mM = 0.7. 

Finally, given our previous simple analytical results of the renormalized self-energies in 
the seesaw limit, see (168]) . (169]) . it is interesting to derive a simple analytical expression for 
the contribution of the heavy neutrino-sneutrino sector to the one-loop radiatively corrected 
Higgs mass in the limit of large mM- Neglecting in (175"]) the contributions from S^- and 
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Figure 14: Contour-lines for the Higgs mass corrections from the neutrino/sneutrino sector 
as a function of the physical Majorana neutrino masses, light \m v \ and heavy m^. The other 
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where Ti V J" denotes the full corrections to the renormalized Higgs-boson self-energy from 



J hh 



the v/v sector and obtained in the mDR scheme as described in the present work. We have 
found that this yields a very good approximation to the full result, i.e. the pole obtained from 
(I75p . In a next step in the above expression E^ has to be replaced by our simplified results 
in the large tum limit, namely, those in (I68bj) and (I69bj) . providing the leading O(m D ) and 
0{m 2 D ) contributions. We have compared numerically this approximate Am™ DR with our 
full numerical results for large trm in Fig. [T2], and found very good agreement, whenever the 
soft SUSY masses are well below %. In fact, the behaviour with vtiu of this approximate 
formula is indistinguishable from the lower line in the left plot of Fig. [12j 
We therefore conclude that the use of the previous ( 1761) with 
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as given in (I68b|) and (I69bj) . respectively, provides an excellent approximation to the full 
result for large Majorana mass values, 10 13 GeV < ttlm < 10 15 GeV and soft masses well 
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Figure 15: Contour-lines for the Higgs mass corrections from the neutrino/sneutrino sector 
as a function of the ratio m^/rriM and the lightest Majorana neutrino mass |m„|. The 
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below rriM, "^susy ~ 10 4 GeV. Furthermore, the above simple approximation can also 
be used for estimates of the differences in the mass correction when applied to the DR 
scheme versus the mDR scheme for different choices of the /zg^ scale. For instance, for 
rriM = = 10 14 GeV and the other parameters set to our reference values as defined in section 

~ l )/M h \ < 1% for 0.1 < fi m /m M < 1. 



14. 2\ we got small differences of |(Am^ 
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5 Conclusions 

In this paper we have presented the one-loop radiative corrections to the renormalized CV- 
even Higgs boson self-energies and to the lightest Higgs boson mass from the one-generation 
neutrino-sneutrino sector within the context of the MSSM-seesaw. The most interesting 
features in this scenario are that the neutrinos, differently to other fermions, are assumed 
to be Majorana particles, and that the origin for the light neutrino mass is not as for the 
other fermions either, but it is instead generated by means of the seesaw mechanism with 
the addition of heavy right handed neutrinos with a large Majorana mass. 

As a first useful result, we have included here the complete set of Feynman rules in 
this MSSM-seesaw context that are relevant for this work, which to our knowledge are not 
available in the literature. These include all vertices for the interactions among the Higgs 
sector and the neutrinos/sneutrinos and for the Z gauge boson and the neutrinos/sneutrinos. 
These Feynman rules have been presented in terms of all the physical masses and mixing 
angles of the particles involved, namely, the CP-even Higgs bosons h and H, the CV-odd 
Higgs boson A, the light and heavy Majorana neutrinos v and N, their SUSY partners v±, 
N± and the neutral gauge boson Z. 

The computation presented here is a full one-loop Feynman diagrammatic one and does 
not make use of any of the approximations applied in the literature. In particular, we do not 
use the mass insertion approximation for any of the involved soft mass parameters, nor we 
neglect the external momentum in the self-energies, which we have found to be relevant for 
the final computation of the Higgs mass corrections. We have presented our analytical results 
in terms of the physical neutrinos, sneutrinos, Z, and Higgs bosons masses. In addition we 
have analyzed the role played by the heavy Majorana mass scale mu, and emphasized the 
differences between the Majorana and Dirac neutrino cases. 

We have fully analyzed the behavior of the neutrino/sneutrino corrections to the renor- 
malized CP-even Higgs self-energies with all the involved masses and parameters: rriM, tan/3, 
Ma, mi, mjj, A v , m v and B v . Our numerical study of the size of these corrections has been 
performed over a wide interval for all these parameters, so that our conclusions can be con- 
sidered as general. From this exhaustive study we have concluded that the most relevant 
parameters are tum, m v, th^ an< ^ Bv m particular, the Majorana mass is by far the most 
crucial one. In general, we have found sizeable corrections to the self-energies, indeed compa- 
rable or even larger than the other relevant one-loop corrections, as the ones from the MSSM 
top-stop sector, at the highest explored values of m M , m u , m^ and B u . We have explained 
here the large size of these corrections in terms of the neutrino Yukawa couplings, which are 
typically large, Y u ~ 0(1) in these seesaw scenarios with heavy Majorana neutrinos. For 
comparison, we have further included the predictions in two renormalization schemes, the 
on-shell and the DR schemes, where we have found interesting differences. These differences 
have been analyzed and explained with the help of simple formulas that are valid in the 
seesaw limit where mu is much larger than all the other mass scales involved. 

The main conclusions on the corrections to the lightest Higgs boson mass are summarized 
in the contour-plots shown in Figs. [H] and [15j For the most conservative scenario of Fig. [Mf 
where all the soft mass parameters are at the TeV scale, the corrections are positive and 
smaller than 0.1 GeV if 10 13 GeV < m M < 10 14 GeV (or, equivalently, 10 13 GeV < m N < 
10 14 GeV) and 0.1 eV < \m„\ < 1 eV. For larger mu and/or \m v \ values the corrections 
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change to negative sign and grow in size with these two masses up to values of around 
—5 GeV for thm — 10 15 GeV and \m v \ = 1 eV. For the less conservative scenario of Fig. [151 
where the soft mass associated to the right handed neutrino sector, m^ is of the order of the 
Majorana mass scale, we find very large negative corrections, at the right upper corner of 
the plot, that is for large tum and m^, of O(10 14 ) GeV, and \m v \ of 0(1) eV. For instance, 
they are around —30 GeV , for m M = 10 14 GeV, m^/rriM = 0.7 and \m„\ = 0.6 eV. In view 
of the anticipated experimental precisions at the LHC and the ILC these corrections are very 
large and should be taken into account if the experimental data indicate the existence of 
Majorana (s)neutrinos. 

In summary, we conclude that the one-loop corrections from heavy Majorana neutrinos 
to the Higgs boson masses are important in this MSSM-seesaw scenario, and overwhelm by 
many orders of magnitude the corresponding corrections in the case of Dirac massive neutri- 
nos. These have also been estimated here and are extremely tiny, smaller than 10~ 22 GeV. 

Finally, we briefly remark on the interesting and more formal issue of decoupling/non- 
decoupling effects from the heavy Majorana neutrinos/sneutrinos sector in the low energy 
MSSM Higgs boson physics. It is clear that our results in the present paper, showing large 
one-loop corrections Am™ DR to the h boson mass for large rriM, suggest that there could 
be indeed non-decoupling effects from the heavy particles in the low energy MSSM Higgs 
bosons physics. Particularly suggesting are the numerical results shown in Figs. 12-15 where 
it is clearly manifested a growing of Am™ DR with tjim- Also our simplified analytical results 
for Am™ DR in (168bj) . ( 169bj) . (|76|) and (177)) suggest a non-decoupling effect, since the mass 



correction does not vanish in the asymptotic limit rriM — > oo, even for Y v (or nip) kept 
fixed. However, we believe that one should not conclude on non-decoupling effects based 
just on the behaviour of the Higgs mass corrections with rriM- It is well known that the mass 
itself is not the proper physical observable to study the decoupling/non-decoupling issue. 
A more proper tool for that study would be the use of Effective Field Theory techniques, 
and more concretely the computation of the one-loop effective action by integration in the 
path integral of the heavy degrees of freedom. An expansion, valid to low external momenta, 
p <C mM, of the derived 1PI renormalized Green functions with Higgs bosons in the external 
legs would provide the definite answer to the issue of decoupling/non-decoupling of the heavy 
vr, z/r, degrees of freedom in the low energy Higgs boson physics. Alternatively one could 
perform one-loop predictions within the present MSSM-seesaw model for other more proper 
observables for this issue like, for instance, cross sections involving Higgs particles in the 
external legs, decay rates of Higgs bosons, etc. The behaviour of these kind of radiative 
corrections at asymptotically large mu could also be conclusive on this issue. All these 
proposed studies are extremely interesting but are far beyond the scope of the present work. 
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Appendix A: New Feynman rules 

In this appendix we collect the Feynman rules within the MSSM-seesaw that are relevant for 
the present work. These correspond to the interactions between the neutrinos and sneutrinos 
with the MSSM Higgs bosons and between the neutrinos and sneutrinos with the Z gauge 
bosons. We write all the Feynman rules here in the physical basis. Here c w = cos 6^. 



h, H, A 




^m D sin20(^,m,-; T5 c O t/3 



h, H, A 




-z^m D sin2£(fff,f^,-z 75 cot/3) 



h, H, A 




-ij^m M sm9cos9 (§- g^,-i 75 cot ft 




^cos 2 ^7m75 




^sin^ 7 , 75 




2^- sin 9 cos 9 7^75 



Table 1: Three-point couplings of two Majorana neutrinos to one MSSM Higgs boson and 
of two Majorana neutrinos to the Z gauge boson. 
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c^ i 4cwM 9 8in ^ -4c ff cosamj + 2c„ cos a m D (A u + m M + /iitana) sin2# H 

u+ X +^f sin/3sin(a + /3)(l + cos2# + )] 

/i 
<-^ —i 4e M g sin a [4c w cos a 77i I, + 2c w cosam£)(y4„ + 777.M + //tana) sin20_| 

JV N s -^sin/3sin(a + /3)(l-cos20+)] 






' v -^ 2cwM ^ sin/3 [c w cos a m D (A, + m M + \i tan a ) cos 26>+ 

s 

N + N — ^sin/3 sin(a + /3) cos 6»+ sin 6»+] 
i?_ • ' 
-c" i 4CwM ^ sin/3 [~ 4c w cos a m l> + 2cw cos a m D (A u - m M + // tan a ) sin 20. 

N 

,}_\ n + ^fsin/3 sin(a + /3)(l + cos2#_)] 

iV_ ,' 

h 
-<' —i 4CvjM 9 w sin a [4c w cos a m 2 D + 2c w cosam D (y4„ - m M + /x tan a ) sin 20. 

]V_ N -^fsin/3 sin(a + /3)(l-cos2^_)] 

h ,''' 

c' -1 2CwM 9 w sin ^ [c w cos a m D {A v - m M + /1 tan a ) cos 26L 

s M 2 



— sin /3 sin(a + /3) cos 9_ sin 0_] 



AL V N 



Table 2: Three-point couplings of two sneutrinos to the Higgs boson h. The corresponding 
couplings to the Higgs boson H are obtained from the ones here by replacing cos a — > 
sin a , sin a — > — cos a , sin(a + /?)—>— cos(a + (3). All the couplings not shown here vanish. 
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/ V- 

-£- cos 9+ cos 6L (p + p% 
p x ^ 



£- sin 9+ sin 0_ (p + p% 



Z, 



iV 4 



^ ^- sin 9+ cos 9- (p + p') i 



P V V 

V- 



Z P ' 
^r^lLj-.s ' 7T- cos 9 + sin 9_(p + p') u 

p ^^ 

Table 3: Three-point couplings of two sneutrinos to the Z gauge boson. All the couplings 
not shown here vanish. 
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N+ N - 



? y 4 



' -< " 2M W 



2/_ 



j^- cot /3 m/) [(Aj, + /i tan /3) sin(#_ — 6 + ) — tum sin(6L + 9 + )] 



Table 4: Three-point couplings of two sneutrinos to the Higgs boson A. All the couplings 
not shown here vanish. 



40 



t 
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h, H #4 

X^ z ' sSvfc?fl [ 4 (~ cos2 a ' sin2 a ) c2 v m D + cos2aM^ sin 2 0(1 + cos 20+)] 

h,H v+ 
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IK -i^^—^[A(cos 2 a,sin 2 a)clm 2 D (-,+)cos2aM^sm 2 /3(1 -cos 29+)] 



u w lvJ -w l 



h, H n + 

h, H v+ 



x x „ (+, — ) i-^pr cos 2a cos#+ sin# H 



ft ^+ 



X V 2 a/ sin 2 8 sin 2a [~ 2 c>d + ^ sin 2 /3 (1 + cos 29+)] 



s >. 



' 8c|M^ sin 2 /3 

H v+ 

ft N+ 



" ' -i ^J sin 2 ^ sm2a[2c 2 w m 2 D - M^sin 2 (1 - cos29+)] 



x' "Se^M^sin 2 /? 



# N+ 
ft K 

>. •■ 

, X N z'/t- sin 2a cos # + sin# + 

# ' * N+ 

Table 5: Four-point couplings of two sneutrinos to two CV-even Higgs bosons. The corre- 
sponding couplings for i>_ and iV_ can be obtained from these by replacing 9+ — > 9-. All 
the couplings not shown here vanish 
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84M|sin^ t 4 cos2 Z 3 c w ™z> - cos 2/3 M 2 , sin 2 /3 (1 - cos 2^ 



iV 4 
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Z„ 
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ij^r cos 2/3 cos# + sin 6^ 



% i~ cos e +9»* 



z ^ sin e +9^ 



2c s 



cos 0+ sin 6+gnv 



Table 6: Four-point couplings of two sneutrinos to two CP-odd Higgs bosons and of two 
sneutrinos to two Z gauge bosons. The corresponding couplings for z/_ and iV_ can be 
obtained from these by replacing 9 + —*•#_. All the couplings not shown here vanish. 
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Appendix B: Majorana case. One- loop neutrino/sneutrino 
corrections to the unrenormalized self-energies and tad- 
poles 

In this Appendix we collect all the analytical results for the neutrino and sneutrino one-loop 
corrections to the Higgs boson tadpoles and unrenormalized self-energies, and to the Z self- 
energies, within the MSSM-seesaw. The contributions from neutrinos iy) and sneutrinos iy) 
are presented separately for clearness. Here c w = cos#w 

+mr,n sin a sin 26 + + mjj cos a (— 2mo + [A v + mj/) sin 29 + )) 
+ A [m?J(Mf cos 2 #_sin/3sin(a + /3) 

+mo/j sin a sin 26L — mo cos a{2rri£, — (A u — tum) sin 28-)) 
- A [m^ ](-M| sin /3sin(a + /3) sin 2 6>+ 



+2mr, cos a(mr> -\ — (A u + m M ) sin 26 + ) + m D fi sin a sin 26 + ) 

- A [m^J(-Mf sin/3sin(a + /3)sin 2 6L 

+2rriD cosa(m£) + -(A„ — mu) sin26L) + mu/xsinasin26'_)) (79) 

^hh(P 2 ) = ~ RA 2 M 2 2 CQS a T R 29 [2m 2 D A [ml} + (2m 2 D + m 2 M )A [m 2 N ] 
64c4iw|7r / sin p 

+ Am 2 D mlB [p 2 , m 2 , m 2 ] + m 2 M {m 2 v + m u m N )B [p 2 , m 2 , m^] 

+ Am 2 D m 2 N B [p 2 , m 2 N , m 2 N ] 

+ /(2m 2 3 fi 1 [p 2 ,m 2 ,m 2 ]+mi f 5 1 [p 2 ,m 2 ,m 2 v ]+2m 2 3 fi 1 [p 2 ,m 2 v ,m 2 v ])] (80) 

E ^(P 2 ) = mo 2 a/2 2 • 2/ , [-4A [m 2 + ](-2m 2 ,cos 2 a + Mlsin 2 /3cos2 a cos 2 g + ) 
51 2c^M^iT z sin p + 

- 4Ao[m^ ](-2m^cos 2 a + Mf sin 2 /3cos2asin 2 6> + ) 

- 4A [m 2 ~_}(-2m 2 D cos 2 a + Mf sin 2 /3 cos 2a cos 2 6L) 

- 4v4 [m| r ](-2m|,cos 2 a + M|sin 2 /3cos2asin 2 6'_)] 

+ 2Bo[p ,m^ ,m s+ J(4m Z) cos 2# + cos a (A v + m^ + /i tan a) 

+M§ sin /3 sin(a + /3) {M% sin (3 sin(a + /3) sin 2 26 + 

—2mr> cos a(v4„ + m M + /x tan a) sin 46 + ) 
+ 2i? [p ,?ti 7 y ; m i>_ J (4m D cos 2#_ cos Q(Aj,-mM + /itana) 

+M| sin /3 sin(a + (5) (Mf sin /3 sin(a + 0) sin 2 20_ 
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—2mr> cosa(74„ — tum + /J tan a) sin46L) 
+ 4_B [p 2 ,m^ ,m 2 ^ ]{m D cosa(2m D + sin29 + (A u + m M + /itana)) 

-M| sin /3 sin(a + /?) sin 2 9 + ) 2 
+ 45 [p 2 ,m 2 y ,771^ ](m D cos a {2m D + sin 29_(A U — m M + /j tan a)) 

-M z sin /3 sin(a + (3) sin 2 6L ) 2 
+ 4:B [p 2 , ml , ml ](mD cosa(—2mo + sm29 + (A u + mu + /itana)) 

-M| sin p sin(a + /3) cos 2 9 + ) 2 
+ 4:B [p 2 ,ml_,ml_](m D cosa(—2m D + sm29_(A u — m M + /xtana)) 

-M|sin/3sin(a + /3)cos 2 ^) 2 ] (81) 

The corresponding results for the tadpole Tjj, and the unrenormalized self-energy T^hh 
are obtained from the above formulas by replacing cos a — > sin a , sin a — > — cos a , sin(a + 
f3) ->■ -cos(a + /3). 

E ^^ 2 ) = - ^2MiJ m2a T R W [2m 2 D A [ml] + (2m 2 D + m 2 M )A [m 2 N ] 
12qc^M z tt z sin p 

+ 4m 2 D mlB [p 2 , m 2 , m 2 ] + m 2 M {m 2 v + m u m N )B [p 2 , m 2 , m^] 

+ 4m 2 ,m 2 v 5 [p 2 ,m 2 v ,m 2 v ]] 

+ /(2m 2 5 5 1 [p 2 ,m 2 ,m 2 ]+m 2 f 5 1 [p 2 ,m 2 , 7^^ + 2772^ [^,771^,771^])] (82) 

S ^(P 2 ) = , 192 J 2 ■ 2 J 4A [m 2 + ]sin2«(m^-M|sin 2 /3cos 2 g + ) 
bl2c^M z 7r z sm (3 + 

+ 4A [m^ ] sin 2a{m 2 D - M 2 Z sin 2 (3 sin 2 0+) 

+ 4 A [ml_ ] sin 2a (m|, - M 2 Z sin 2 /3 cos 2 0_ ) 

+ AA [m\_ } sin 2a (m 2 ^ - M| sin 2 (3 sin 2 #_ ) 

+ 2B [p 2 ,m% + ,ml + ] x 

(2m|) cos 2 29 + (—2(A u + tum)^ cos 2a + ((A^ + mM) 2 — l^ 2 ) sin 2a) 
+M| sin/3(-M| sin/3 sin(a + /3) cos(a + 0) sin 2 26> + 
+m D ((A„ + m M ) cos(2a + /3) + /isin(2a + (3)) sin46» + )) 

+ 2B [p 2 , m%_,ml J x 

(2ml, cos 2 29_(—2(A u — m A ,/)/i cos 2a + ((v4„ — m M ) 2 — fi 2 ) sin 2a) 
+M§ sin /3{-M z sin /3 sin(a + (3) cos(a + (3) sin 2 20_ 
+m/j((74 ;/ — itim) cos(2a + (3) + /isin(2a + /?)) sin46L)) 

+ 2B [p 2 ,m 2 ^ ,m 2 jy }(rn 2 D (— 2ji cos 2a sin 29 + {2m D + (A u + m M ) sin 26' + ) 

+ sin2a(4m|j + 4m D (A v + m M ) sin 2^+ + ((A„ + m M ) 2 - /J 2 ) sin 2 20+)) 
+M z m D sin /3 sin 2 0+ (2/x sin(2a + /3) sin 20+ 
+2(2m D + (A, + m M ) sin 20+) cos(2a + (3)) 
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Xaa(M 2 a ) 



-M z sin 2 sin 4 0+ sin 2 (a + 0) ) 
+ 25 [p 2 ,m 2 y ,m^v ](— m^^/i cos 2a sin 20_(2m D + (A u — m M ) sin20_) 

- sin2a(4m|, + Am D (A u -m M )sin20_ + ((A„ - m A/ ) 2 - /i 2 ) sin 2 20_)) 
+M|m D sin /3 sin 2 0_ (2/i sin(2a + /3) sin 20_ 

+2(2m D + (A„ - m M ) sin 20_) cos(2a + /?)) 
-M| sin 2 /3 sin 4 0_ sin 2 (a + /?)) 
+ 2i? [p 2 , , m. 2 , m 2 + ](— m 2 : )(2/icos2asin26' + (— 2m/) + (A/ 

- sin 2a(4m 2 D - Am D (A u + m M ) sin 20+ + ((A u + m M ) 



ibm) sin 20+) 
/i 2 ) sin 2 20+)) 



+M 2 z m D sin /3 cos 2 0+ (-2/i sin(2a + /3) sin 20^ 



+2(2m Z) - (A, + m M ) sin 20+) cos(2a + /?)) 
-Af| sin 2 /3 cos 4 0+ sin 2 (a + /3)) 
+ 2i?o[p 2 ,m 2 _,m 2 _](— m|,(2/icos2asin20„(— 2mz) + (A v — ttim) sin20_ 



— sin2a(Am 2 D — Amr>{A u — tum) sin20_ + ((A 
+M 2 z m D sin /3 cos 2 0_ (-2/i sin(2a + /3) sin 20_ 
+2(2m D - (A - m M ) sin20„) cos(2a + /3)) 
-M z sin 2 /3 cos 4 0_ sin 2(a + /3))] 



r«M 



/i 2 )sin 2 20_)) 



cos 2 /3 sin 2 20 



64cIM z tt 2 



sin /3 



[2m^A[m 2 ] + (2m 2 D + m 2 M )A [m 



N\ 



13) 
34) 



m 2 M (ml - m u m N )B [M\, m 2 , m^] 
M 2 (2m 2 9 (E 1 [M 2 ,m 2 ,m 2 ]+5 1 [M 2 ,m 2 v ,m 2 v ]) + m 2 /J B 1 [M 2 ,m 2 ,m 2 v ])] 



S^(M 2 ) 



2564M|tt 2 sin 2 /3 



[A [m? ] (4m|, cos 2 /3 - 2M| cos 2/3 sin 2 cos 2 0^ 



+ A [m% } (Am 2 D cos 2 /3 - 2M| cos 2/3 sin 2 (3 sin 2 H 
+ A [m~_ ] (4m^ cos 2 /3 - 2M| cos 2/3 sin 2 cos 2 0_ 
+ A[™i }{Am 2 D cos 2 /3 - 2M| cos 2/3 sin 2 /3 sin 2 0^ 



+ 4m 2 ^[ J Bo[M 2 ,m? + ,m?_](/isin/3sin(0_ 
+ cos /3( A sin(0„ — 0+) — m M sin(0_ 

+ Bo[M%m\ ,m 2 ? _]( ) wsin/3sin(0_ -9. 
+ cos /3( A sin(0_ — 0+) + m M sin(0_ 

+ B [M a, rn 2 ^,ml + }(fi sin /3cos(6^ — 4 

+ cos /3( A cos(0_ — 0+) — rriM cos(0_ 
+ B [Ml, m% + ,ml_](Li sin j3cos(6_ - 4 

+ cos/3(A cos(0_ — 0+) + m M cos(0_ 



^ + ))) 2 
+ ))f 
+ ))f 

o + ))Y]\ 



55) 



£zz(Ml) 



32ci7r2 



cos 4 A [m 2 ] + - (3 + cos 20) sin 2 9A [m 



N\ 



45 



M 2 
+ 2cos 4 # ( mlB [M*,ml,ml] - B Q0 [M 2 ,ml,m 2 v \ + -^B x \M 2 z ,m 2 v ,m 2 v _ 



( M 2 

+ 2 sin 4 9 \m%Bo[M%,m 2 N ,m 2 N ]- Bqo[MI, m 2 N ,m 2 N ] ^ -^B x \M 2 z ,m 2 N ,m 2 N \ 

+ - sin 2 19 (m u (m u + m N )B [M 2 z , m 2 , m 2 N ] - 2B 00 [M Z , m 2 , m 2 N ] 

+M 2 z B l \M 2 z ,m 2 u ,m 2 N })] (86) 

X ZZ (M 2 ) = ^L- 2 \A Q [m 2 J cos 2 9„ + A [m 2 + ] cos 2 9 + 

+ A [m^_ ] sin 2 9- + A [m\ + } sin 2 9+ 

A(B 00 [M 2 Z , ml + ,ml_] cos 2 9_ cos 2 0+ + B 00 [M 2 Z , m%_,ml + ] cos 2 #+ sin 2 9_ 

+ S o[Ml, m^ + , m 2 _] cos 2 0_ sin 2 # + + 5 00 [M|, m^ + , mjj sin 2 #_ sin 2 0+)] (87) 

The definitions of the loop functions Aq, Bo, B\ and -Boo appearing in this and the next 
appendices can be found, for instance, in Ref. [IS] (where Bq§ = B22). 

Appendix C: Dirac case. One-loop contributions from 
neutrinos and sneutrinos to the renormalized h Higgs 
boson self-energy 

We present here the result for the one- loop corrections from neutrinos (u) and sneutrinos 
(i>) to the renormalized hh self-energy in the case of Dirac neutrinos, obtained in the DR 
scheme. Here c w = cos 9w- 

SL^Wc = 32/^2^2 {^oK](sin 2 (a + /3)M 2 

'sin(2a - 3/3) + 3 sin(2a - B) - 2 sin /3))m 2 D 



sin/3 
+ sin 2 (a + B)M z {m 2 D B [M z , m 2 D , m 2 D ] 

-2£? 00 [M 2 , m 2 D , m 2 D ] + M 2 Z B X \M 2 Z , m 2 D , m 2 D }) 

2 

C-t , a r O O O -\ O r ~\\ 

- 2 . 2 n ( 2m D B o[P ,m D ,m D \ +p Bi\p ,m D ,m D \) 
sin p 

2 , cos 2 (q - /3) cos 2 /3 2 , , -> 

+ 2M A m D r^— #1 M A , m D , m D f 

sin o J 



sin (3 



,2 



. r 9 1 r , ,9 9/ ^ 9^ sin(a — 8) sin 20 
E^(p 2 ) Dirac - — l— 2 J Aolm 2 ] 8M| sm 2 (a + B) cos 2 + 2m D l , ^ x 



-w J "Z' 



sin/3 



(/i(3sina — sin(a — 2/3)) + A v (3cosa + cos(a — 2/3))) 
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+ A [ml 2 ] 8M| sin 2 (a + (3) sin 2 9- 2m D 



sin(a — (3) sin26 l 
sin/3 



x 



(/i(3sina — sin(a — 2/3)) + ^4^(3 cos a + cos(a — 2(3))) 
1 1 



;S [p ,m iil ,m c 



16sin 2 /3~ L 

r2 



2(8m 2 3 -M|)cosa; 



2M|(cos(a + 2(3) - 2 cos 2# sin /3 sin(a + (3)) 



-8rriD sin 26* cos a (A u + /^tana) 



lb sin p 

r2 



2(8m|, -M|) cos a 



1 



8 sin 2 (3 



2M|(cos(a + 2/3) + 2 cos 20 sin (3 sin(a + /3)) 



+8rri£, sin 2$ cos a(A^ + fi tan a) 



B [p 2 ,ml 2 , m 2 J - 4M| sin 29 sin /3 sin(a + (3) 



+8rriD cos 26 1 cos a(A v + /x tan a) 
+ 8m 2 } cos 2 (a-/3) cot 2 /3S [M|,m 2 2 ,mlJ (A, + /i tan/3) 2 

8A/f sin 2 (a + /3) (2 cos 4 ^ 00 [M|, m 2 m 2 J + 2 sin 4 ~BB m [M 2 Z: m%^m\ 



+ sin 2 2QBw[Ml,m\,m\ 



39) 
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